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PREFACE 


Having served a number of years in the capacity of communications 
officer in U.S. Marine Corps Infantry Divisions, the writer has been 
aware of the need for a method of analysing military communication 
systems. The delay of messages has been and still is of major 
concern to all military communicators. The increasing speed of action 
of military units ages information rapidly; and as a result message 
delays achieve greater significance. 

It was thus with great pleasure and interest that the writer 
undertook the problem of investigating military message delay times 
as recommended by Dr, Edgar Reich of The Rand Corporation. The major 
part of the work was accomplished at The Rand Corporation, Santa 
Monica, California during the industrial tour phase of the Electronics 
Curriculum as presented at the U.S. Naval Postgraduate School, 
Monterey, California. 

The writer takes pleasure in acknowledging his indebtedness 
to Dr. Edgar Reich, whose assistance, encouragement and instruction 
in queueing theory were prerequisites for the preparation of this 


work. 
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INTRODUC TION 


The problem of formulating a systematic method for the determina- 
tion of communication requirements of Army combat units is being 
investigated by Haller, Raymond and Brown, Inc. E In determining 
requirements, it is first necessary to establish relationships 
between the communication load and message delay times for various 
communication systems. Commercial communication organizations have 
in the past applied queueing theory to the solution of similar problems 
fe However, since the communication load and discipline of the 
military are different from those of civilian users, the presently 
available tools of queueing theory are not directly applicable to 
military problems. Research in two directions is thus required. It 
is necessary to determine more accurately the load characteristics of 
military communications and to further develop queueing theory with 
military applications in mind. This paper discusses some of the 
queueing theory techniques available and develops the theory a step 
further, 

Queueing theory is the application of the mathematical theory of 
probability to problems involving traffic congestion i31. Tha traffic 
congestion occurs when a service facility is subjected to service 
demands which vary statistically. The traffic may be aircraft attempt- 
ing to land at a airport, customers at a bank teller's window, or 
messages in a communication system. ‘Ve assume that there ara sig- 


nificant statistical fluctuations of either the inter-arrival times, 





or the service times, or both; and that it is undesirable for the 
traffic to wait. An analysis of the waiting lines or queues that 
result may lead to methods for shortening the lengths of the lines and 
thereby decreasing the delays encountered by the traffic. 

Chapters II and III contain the necessary groundwork for an under- 
standing of the application of queueing theory to communication 
problems. | 

Chapter IV contains the solution to a classical queueing problem. 

Chapter V contains the major contribution of this work, consist- 
ing of the solution to the queueing problem characterized by the follow- 
ing: 

a. single channel 

b, messages from an exponential length distribution 

Ce messages arriving in accordance with Poisson Laws 

d. messages transmitted in accordance with assigned priorities. 

The solution consists of mathematical expressions and graphs 
from which the communicator is able to compute the probabilities that 
messages of the various priorities will be subjected to particular 
delays. It is thus possible for him to determine how well the channel 
will perform and perhaps make any necessary changes to correct un- 


desirable delays. 








CHAPTER II 


PARAMETERS OF THE COMMUNICATIONS CHANNEL 


ime Introduction 
The over-all problem of military communications involves the ex- 
change of information among a number of sources. The amount of ine 
formation and tne number of sources depend on the size of the unit 
and the tactical situation. The solution of the problem would result 
in determining the number of channels to be established and communica- 
tion procedures to be enforced. It appears then, that the solution 
to the over-all problem would require an extensive knowledge of the 
influence of various parameters on the operation of a single channel. 
Single channel as used here-in applies to a standard military radio 
net. 
In order to treat the operation of a channel mathematically, it 
is necessary to specify the following: 
l. Message generation times 
2. Nessage length distribution 
3. Zueue discipline 
As noted earlier, either the message generation times or the 
message lengths, or both, fluctuate statistically. The application of 
queueing theory may then result in the calculation of the probability 
that a message will wait a time t between it's origination ty the 
sender and the beginning of it's transmission over the channel. This 
probability is termed the delay or waiting time probability. It should 


be noted that the time t does not include the transmission time of the 





message but is the time the message waits until it's transmission begins. 
2. Message Generation Times 

Message generation times refer to the intervals of time between 
the origination of messages by the units utilizing the channel. These 
time intervals may be described by stating the probabilities of having 
various numbers of messages generated per unit of time; or the proba- 
bilities that intervals between message arrivals will be of various 
lengths. The message generation times for military units will vary 
depending on the size of the unit, the mission, and numerous other 
factors. Haller, Raymond and Brown, Inc. is attempting to determine 
empirically the message generation distribution for company-size units. 
There &re undoubtedly dependencies between messages since incoming 
massages prompt outgoing messages. However, in large units it is 
expected that these dependencies average out. 
Se Message Lengths 

Message longths may be described by stating the probability that 
a message chosen at random will have a length less than a given amount. 
The message length distribution may also be influenced by the tactical 
situation. 
4. Queue Discipline 

Tha term queus discipline refers to the handling of messages 
awaiting transmission, Military communications procedure usually 
requires that delayed calls be serviced consecutively in the order of 


their origination unless a priority system is in effect. 





CHAPTER III 


PROBABILITY DISTRIBUTIONS OF THE PARAMETERS 


introduction 

The exact probability distributions of msssage generation times 
and message lengths of military units have not as yet been deternined. 
However, if the discussion is limited to units of about division size 
during, periods of greatest activity, it may be assumed that the message 
length distribution is exponential and the message generation times 
distribution (number of messages generated per unit time) is Poisson. 
These are considered the vest assumptions until empirical data de- 
monstrates otherwise. 
2. Exponential Message Lengths 

Suppose that & message is being transnitted, then message length 
and transmission time may be used synonomously. Let the length of 
transmission , in seconds, be represented by a random variable t; 
possessing a prohability distribution. Th? transmission may end during 
any small interval of time. If the probability of termination during 
any small time interval is constant and independent of how long the 
transmission has been in progress, the message length distribution 
is exponential. Telephone company data demonstrates that the ex- 
ponential distribution is an accurate portrayal of telephone convarsa- 
tion lengths [4]. 

The exponential distribution may be derived from the assumption 


that the probability of & message having a length between t and t +h 


>. 





(h is a short tine interval) is constant and independent of t (4) L5) . 
Here we will accept this distribution function and investigate its 
characteristics, 
If a message is chosen at random fron an exponential distribu- 
tion, the probability that its length is less than t is 
nase Icy 


The message length probability density function is 


ORO ND DE Eo 
dt 


The average message length is 


CK OO 
le lee tuerAtat = 1 
5 JE 


Let zero time represent the time at which transmission of a 
randomly chosen message begins. Let t represent a time of observation 
at or subsequent to time zero. Then to denonstrats that the probability 
that t!e transmission will end during h is independent of t, write the 
probability that a message be of length between t and t +h 

Ner EIE) 
Let P = conditional probability that if a transmission has been in 
progress for a time t, the message length is between t and t +h. 


Then by Bayes! Theorem 
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Thus if a transmission has been in progress for a time t, the probability 





that it will end during the following interval of time h is 
ie we 
and is independent of t. 
If h is assumed to be small, then by expanding the exponential and 
discarding terms in h* and higher powers 
LE een 
5. Poisson Message Generation Times 
If messages are generated independently of each other, and if 
the probability that a message is gensrated in a time interval is in- 
dependent of the number of messages generated in previous time intervals, 
the message generation process is considered to be Poisson. The Pois- 
son distribution may be derived by postulating that the probability of 
an event occuring in a time interval is a constant and independent of 
the number of events that occurred previously [4] si. Here we accept 
this function and investigate its properties. 
For Poisson arrivals the probability that n arrivals occur during 


a time interval t is 


r.(t) z (A)! e ^t 


The average number of arrivals per unit time is 


ey ale) „aA re 
ic ar tn} 

nzo nzo 
If h represents & small time interval, then 


fo(h) s€7 Anzi -Àn 


is the provability that there will be no arrivals during time h; and 





fj (h) = Ah eo ae 
is the probability that there will be one arrival during h. All pro- 
babilities of more than i arrival consist of terms of order h* and 


higher and are discarded. 3 





CHAPTER IV 
EXPONZNTIAL MESSAGE LENGTHS, POISSON ARRIVALS 


AND ORDERED QUEUEING IN A SINGLE CHANNEL 


me introduction 

Assume that the population of messages to be transmitted pos- 
sesses an exponential length distribution, the message generation 
process is Poisson, and that the messages are serviced ona "first 
come, first served" basis. The information desired is the distribu- 
tion function of the waiting time. The attack on the problem consists 
of determining the probability that the channel is in a particular 
state and the conditional probability that a message will wait less 
than time t if arriving while the channel is in a particular state, 
The summation of the products of these two probabilities will result 
in the probability that a message will wait a time less than t before 
its service bagins. 
2. State Probabilitiss 

Define system state n to be the state of the channel when there 
are n messages in the system being served or awaiting service, P) 
is used to designate the probability that at a time t the system state 
is n. If h represents & small interval of time, then from the discus- 
sion on distributions, wh represents the probability that & transmis- 
sion will terninate during h and Ah represents the probability that 
a message will arrive during h. 

List now the possible ways in which the system can arrive at state 


n,(n>1), at time t: 





Im ac Ces hn the state isn., no arrivals and no transmission 
completions during he 
2, at t - h the state is ntl, no arrivals and one transmission 
Completion during he 
oe at t - h the state is n - 1, one arrival and no trans- 
mission completions during he 
4, all others (states other tian n, n-l, n4l, and multiple 
arrivals and /or completions). 
The probabilities of these events are: 
1. (1- 45) 1- AR), (6-28) z- (o 4 A)n] P (t-h) 
2. (4 h)(1- Ah)Pa qu (t-h) =unP, , y (t-h) 
3. (Ah)(1-ph)P,_) (teh) % AP, _, (t-h) 
4. order of h* and higher 
The sum of these is the probability that state n will be reached 
at time t, 
PA (£)e [1-Gt 9A )h ] P4 (&-h) e p hP, (62h) AnP, 4 (t-h) 


p t)-P (t-h)= = +pL)P (toh)+pP (toh)+ AP (t-h) 


-l 





h 
lettinag h0 


(1) ap EEA AE E pr, (RAR (t) mA 


-1 





at 
Now consider state zero: 
1. at teh state is zero, no arrivals and no transmissions 
completed during h 
2. at t=h state is l, no arrivals and one transmission com- 


pleted during h 


LO 





Se all others 

Incenrobabilities are: 
l. (L-uah)(1-A h)P (th) = |1-(z + A Ja! PT 
cane h)(1- Ah)P, (t-h) «jh, (t-h) 


2. order of h* and higher 
Po(t)s [R-i & A )5] Pg (7h) phP, (th) 
Po(t)-Po(teh)= -(u +A )P (t-h)+ uP, (th) 


h 


letting h70 


(2) ap (t)= -( € A )P Ct) eu P, Ct) 
dt 





It has been demonstrated [6] that when À< u and t is large, that 


the effect of the initial conditions are overcome and that the 


derivatives of equations 1 and 2 approach zero. Then 
(A +p)P = pe tar 0M 
AP SH D 
Assume a solution of the form P six. Then 


(A $42) ¥ =x 4) 


xs A ,x=l 
co y 
Since ) P (t)=1 for all t, the solution x#l is discarded, 


n=O 
Le 








and xa A, then 
44 i 
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n 
A 
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3. Waiting Time Distribution 

Let W (t) represent the conditional probability that if a message 
arrives when tne channel is in state n, its waiting time is less than 
ips c bet w(t) be the probability density function. Consider a message 
M arriving when the channel state is n. There are then n messages 
each with an exponential length distribution to be transmitted prior 
to M, The density function of M's waiting time is the convolution 
of n exponential density functions 

ne DD 
nel). 
The density function of a message's waiting time when the 


channel state is unknown is , 


DO 


w(t)z y w (t) ME" "e à 2 uer LOT 
Gere 
na0 
CKD 
sett. A) (At)? =AQ- buon 
fe ^ni 
n=0 


12 





The distribution function of tho waiting time is 


W(t) fw(t)at 


E D» fe-(u -À tat 


`~ 


Applying the condition that "(O)z(l- =). since the probability 
of waiting zero time is merely the probability that the message 


arrives when the channel is in state zero, 


Y(t)=1- a ar he 


The mean waiting time may be computed from 


I 
t- few(t)at 
O 


E A) cea At 


O 


t= À 
Ag za) 
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CHAPTER V 
POISSON ARRIVALS, EXPONENTIAL MESSAGS LENGTHS, 


AND PRIORITY QUEUEING Iii A SINGL® CHANNEL 


hee lntroduc tion 
The problem solved in Chapter IV is among the simplest in 
queueing theory 17, Among the problems that have been solved are: 

Single Channel 

l. Poisson arrivals, exponential serving tine, service 
in order of arrival. 

. 4,Poisson arrivals, exponential serving time, random 
service. 

ó, Poisson arrivals, uniform serving time, service in 
order of arrival. 

4. Poisson arrivals, Pearson 111 serving time, service 
in oọorsr of arrival, 

o, Arrival rate dependent on channel state, Pearson III 
serving time, service in order of arrival. 

6, Poisson arrivals, any serving time distribution, 
service in order of priorities, solved for mean wait- 
ing time only. 

Multiple Channel 

l. Poisson arrival, exponential serving time, service in 
order of arrival, 

2. Poisson arrivals, exponential serving time, service in 


order of priorities, solved for mean waiting time only. 


14 








The two priority cases mentioned above are of interest to 
military connumicators. However the solutions exist for mean wnit- 
ing tine only and furnish merely an indication of the delays encounter- 
ed by the messages. In military operations messages subjected to 
very long delays &re not compensated by messages having very short 
delays, The military communicator would prefer knowladge of the 
waiting time distribution function which would enable him to evaluate 
his communication system complotely. In this chapter a waiting time 
distribution function is derived for a single channel characterized 
by the following: 

l. Poisson arrivals 

2, exponential serving time 

Oe Service in order of priorities 
2e Assumptions 

The following discussion is limited to the case in which the 
messages arrive independently of each other and the arrival proba- 
bilities are independent of the state of the channel. Thus the units 
arrive in accordance with independent Poisson Laws with Ap represent- 
ing the average arrival rate of priority p. It is also assumed that 
all priorities possess identical exponential message length proba 
bility distribution functions. 

The number of different priorities is Y , with l representing 
the highest and r the lowest priorities, The queue discipline follows 
standard military communications procedure which results in higher 
priority messages displacing messages of lower priority in the waiting 


line (but no disturbance of the message being transmitted), and 


15 





messages of equal priority being transmitt-d in order of arrival. 
Imwprineinle the attack on the problem consists of: 
1. deternining the conditional density function w (t) 


n,p 


2. determining P 
n,p 


3. determining the unconditional density function o by 


co 
) or w (t) 

n,P n,p 

n=0 

$. Conditional "aiting Time Distribution 

In attempting to deternine the conditional probability t'at a 
message of priority p will wait t or less if arriving when the 
channel state is n,p, consider a message Y of priority p arriving 
when the system is in l,p,i.e., no messages of priority p or higher 
awaiting service and one message of any vriority being transmitted. 
During the latter message's transmission, messages of priorities 
higher than p may arrive; and these will be served before X, M will 
wait a timev when the transmission of a message will just have been 
completed and there are no messages of priority higher than p await- 
ing service. The density function of T ("busy period") has been 
determined by Kendall 18] and Bailey [2] , and when applied to the 


present problem becomes 


mac AAA 
‚pP A s EST ^r 


=, 


ssl is used here since the only messages that affect M's waiting 
E 


time are those of priorities higher than p. Since the individual 


16 





arrival distributions of these messages are independent Poisson, the 
sun is also Poisson with mean A s 
2p-l 
Now asstme that a message M of priority p arrives when the 
channel state is n,p. M must wait a time E until the following 
mossages are transmitted: 
l. the n messages contributing to state n,p when M arrived. 
2. 81l messages of priorities higher than p arriving before 
M's service begins. 
Consider an identical system in which a message V! of priority 
p arrives when the channel state is (n-1), p. Y and M’ arrive at 
the same time. Let M' wait t then M must wait an additional 


nel? 


tino E. Thus if d is M's waiting time: 


t at + 
n n-l 


A 


aor for nz2, this recursion tormes 


z - “y 
tzt tr but ENS e 


then the density function of t, : 


SE mo 


and in general, the density function of the waiting time for a message 


of priority p arriving when the channel state is n,p 


wo (tw (t)™ 
rae 1,p 
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Transform pair 576.1, Campbell and Foster uio) furnishes the 


characteristic function of Ww, E 
3 


4 Lu. 


E "Lh a 
(ser y]? 
Where NE Doc E 


- p UESL S 
e Asp FEN Aue P 


The characteristic function of w ‘Ges 
n,p 





C = an” 
n,p l L 
E Hs)? |” 


and utilizing pair 576.1 again 


n n Apr )t I (2t \ À pi) n7o 
(t) C La | "er 


D. " 


) 
n,p E 


0, (6) $ (8-0) 


4. Channel State Probabilities 


In determining the channel state probabilities, note that the 


probability of gtate zero is 


18 





State 1,p can exist in the following ways: 
l, only one message in the channel of any priority 
2. more than one message in the channel all of lower 
priority than p 
ó, more than one message in the channel, the one being 
served of priority p or higher and all others lower 
than p 
The probability of 1,i.e. of the channel containing only one 


message is 


(i= A „(A ) 


"m 


The following may occur resulting in the existence of state 
1,p at time t: 
l. State at time teh is 1,p, no completions, no arrivals 
of priorities p or higher, and either arrivals or no 
arrivals of prioritiss lower than p during h. The 


probability of this event is 


FLA a Py (teh) 


2. State at time t-h is l,p (less the condition of exactly 
one message in the channel), and a completion during 
h. The probability of this event is 


j E 


N A ne 
un|P, (th)-(1- a E 


/ ti 


3e State at time t-h is zero and an arrival of any priority 


13 


during he The probability of this event is 
AnP (t-n)a(l- A Jin 
O m 


4. State at time t-h is 2,p and a completion during h. 


The probability of this event is 


Eu tcn) 
Thus: 
y A n ESSA 
Py p(t ar, Ci) 1 TE AS n)-(1 exp 
Al. , 
4 Anji- me tt h) 


B (t)- P, _(t-h) 


1 
oe thru. (teh 
- AP. oe A ) 
as h0 
A 
— Asp? pit) + MP2 p(t>) 


after initial conditions effects overcome: 


E Seite 
Ast yp? Fe, 
For n?l, the steady state equation is the same as for tho earlier 


case of exponential message length distribution and Poisson arrivals 


except that here for priority p substitute A rA 
2P 


(A +u)P zu? +A P ny! 
s P np (n+1),p ¿p (n-1),p 


Assume as a general solution for n>l 
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Then since TR 


= 
nyo 
IL EK 


P=(1- —) 
O 


Multiplying P and C and summing to obtain the 
n,p n,p 
characteristic function of the waiting time 


mE = AA 
Cal] Ay, A Te am 
p M AL Asp Load st )«4(seq) A 


nal 


A 24 E. | y 
0 A r Mr ter A, ee d, Ge 
p 


: En 
TH Z2 mS 
Á t4 *z NA) "(st A ru -2\ J ju 
2p»-i ` $p-l ip-l' "sp-1l 


At this point it is convenient to convert time to units of 


eL 





average message length and to let 
f A «4 A: 44 
Zp =l sp 
then Azp A =a f "NS 
gp-l £p 


The characteristic function then becomes 


np) 4 —2°U-2F) 
d i E 
(staP Q-2AP )«(sexf ra lsr afalar )? 


The inversion to the density function-can be accomplished but results 


in an infinite sum of Bessel functions 


AORTA f 42) E a E po nlat var ) 
e nal 
However the density function of the waiting time of the highest 


priority messages can be obtained by setting 20. Then 


Fü-^A4P) 
=(1- P)+ (s+1-2P ) 


(021-0 Paare SF DF 


The distribution function of the waiting time of the highest priority 


messages 
t 


"C -1_ pa ~\1-4F jt 
w (t)z as e 
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Figures 1, 2 and 3 are plots of the waiting time distribution 
of highest priority messages for various values of / andÁ. 

An approximate solution a may be obtained by deriving 
the central moments of the waiting time distribution function from 


C,, and utilizing the Edgeworth Series[11 [12], 


y 
l ü-4P)ü-«P) 


Up m 2 2 3 
a-6p) ü-«pP)'  (Q-4P)0-«P) 





>) a Ls RAP) | 
5 3 3 2 4 | 5 
üQ-4P)'"üa-«P) ü-ée)*a-«p)* Q-2P)0-xf) 
- ee 4 AXP (24307 ) O 
4 4 4 3 5 2 6 
au-5P) ü-«P)  üü-56P) ü-«P)  Q-6P) ü-«P) 
; «f 51e af ee: ia 2 
Q- 4P)0-x2)' 
i Hs 1 Ma 
oe Bee (3) 
v (S)zGUz 960- x PA qo 
10 f 
+ o 63 9 Gua. p) 
61 43 
Where 
x=t- Y, 
= 
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Le 


O (x)= Tm e 2 


d qx) 
q^: TE 


AE or 


"n os ),-7; 
A 3 EE. 
BS Vs) eu 1 


Py e v 


l 


The effects of variations of « ,/óand P on the average waiting 
time and the standard deviation of the waiting time are demonstrated 


in Figures 4 and 5. 
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CHAPTER VI 


CONCLUSIONS 


1. Summary of results 


In Chapter V expressions for the following were derived: 


These 


Be 


b, 


d, 


probabilities of tho system states 

the characteristic function of the waiting time 
distribution function 

the probability density function of the waiting time 
distribution function 

the waiting time distribution function of the highest 
priority messages 

an approximate solution for the waiting time distribu- 


tion function of messages of any priority, 


solutions are based on the following assumptions: 


a. messages of each priority arrive in accordance with 


b, 


Ce 


d. 


independent Poisson Laws 

all priorities possess the same exponential message 
length distribution 

the system consists of a single communication channel 


standard military communication procedure is in effect. 


2. Applications o 


In the light of present knowledge of the subject, the assumptions 


of Poisson arrivals and exponential message length distributions are 


considered satisfactory. The assumption concerning the equality of 


the mean message lengths of all priorities is not exact since in 
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general the higher the priority the shorter the message. However, 
if the lowest priority message length distribution is assumed for 

a priorities, the computed delays will be greater then those 
actually encountered and the system should perform at least as well 
2 anticipated. 

The military situation dictates the maximum message delay 
tolerable for each priority. The derived expressions enable a 
military communicator to analize a communication channel and determine 
if it satisfies the tactical requirements of the moment. For example 
the situation may require that at least 99% of the highest priority 
messages be delayed no longer than a particular time. Assuming that 
the assumptions are valid for the case considered, the communicator 
is able to determine how well the channel meets this requirement. 
Increasing the communication facilities, modifying the priority 
system, or further training of communications personnel (to decrease 
average transmission tine) are indicated when requirements exceed 
channel capabilities. 

Military communications is only one application of the results 
presented in Chapter V. Among others are aircraft landings, treat- 
ment of patients, and maintenance facilities. The various priorities 
may arise as a result of emergencies. Needless to say, it is 
necessary that the basic assumptions concerning the distributions 
of the arriving units be satisfied. 

ö. Further Investigation 
A part of the original material presented in this paper resulted 


from an extension of an approach to the priority problem taken by 
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Alan Cobham ¡13] . It is hoped that this paper will be of assist- 
ance to future investigators who attempt to solve the more general 
problems in which the units of the various priorities possess 


different length distribution functions and the number of channels 


is greater than one. 
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